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COSMIC RAY HYSTERESIS AS EVIDENCE FOR TIME-DEPENDENT 
DIFFUSIVE PROCESSES IN THE LONG TERM SOLAR MODULATION 

J.J. O' Gallagher 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland (USA) 20742 

A simple one-dimensional time-dependent diffusion-convection 
model for the modulation of cosmic rays is presented. This 
model predicts that the observed intensity U(t 0 ) at a given 
time t Q is approximately equal to the intensity given by the 
time independent diffusion convection solution under inter- 
planetary conditions which existed a time x in the past, 

[U(t Q ) = U s (t -t) ], where x is the average time spent by a 
particle inside the modulating cavity. Delay times in excess 
of several hundred days are possible with reasonable modulation 
parameters. Interpretation of phase lags observed during the 
1969-70 solar maximum in terms of this model suggests that the 
modulating region is probably not less than 10 a.u. and maybe 
as much as 35 a.u. in extent. 

1. Introduction . Observations of the time behavior of modulated cos- 

mic rays have revealed that the integral intensities above different energy 
thresholds are not modulated in phase, particularly near times of solar 
maximum or minimum (Simpson 1963; Simpson and Wang, 1970; Balasubramanyan 
et al. , 1968; Kane and Winkler, 1969). That is, regression plots between 
the modulated intensities at different energies sometimes exhibit a 
double-valued "loop" which recently has come to be called the "hysteresis 
effect" in cosmic rays, A number of recent studies of this behavior during 
the past solar maximum in 1969-1970 have shown that the effect exists for 
electrons as well as nuclei (Schmidt, 1972; Burger and Schwannenberg, 1973), 
and that protons and helium nuclei are affected differently in a way which 
depends on particle rigidity (Van Hollebeke, Wang and McDonald, 1972 and 
Rygg, O'Gallagher and Earl, 1973). 

Some success in explaining this behavior in terms of the diffusion 
convection modulation model (Parker, 1963, 1965) has been achieved by 
invoking a diffusion coefficient k whose dependence on particle magnetic 

rigidity R varies with time (Schmidt, 1972) or with heliocentric radius 
(Burger and Schwannenberg, 1973) or both (Van Hollebeke, et al. , 1972). 
However as has been pointed out by Rygg , et al. , 1973 all such models are 
fundamentally equivalent providing no observational characteristics which 
can distinguish between them. Furthermore, all involve the introduction 
of several new free parameters in order to provide sufficient flexibility 
to fit the wide range of observed behavior. It is fundamental to each of 
these approaches that the observed spectrum at a particular time is the 
solution to a Fokker-Planck equation under equilibrium conditions i.e. a 
stationary solution in which all time dependent effects are neglected. 

However since the periods of solar maximum and solar minimum are 
probably not equilibrium situations it is instructive to consider an 


1 



2 

177 


analysis based on an alternate model in which the effects of time dependent 
diffusion are included as a first order correction. The model is very sim- 
ple, being based on one dimensional time dependent diffusion which neglects 
all energy loss effects. However the basic conceptions are easily visualized 
and lead to the interpretation of "hysteresis" as a real rigidity dependent 
time lag in the modulated intensity at different rigidities. The model in- 
volves no new parameters and provides a basis for using hysteresis to mea- 
sure some of the dynamical characteristics of the modulating region. In 
particular, when applied in a simple minded fashion to the balloon observa- 
tions of hysteresis reported by Rygg et al. , (1973) it provides a crude 
estimate of the distance in the boundary of the modulating region of between 
10 and 20 astronomical units (a.u.). However when the limiting behavior at 
small diffusion coefficients is taken into account the model suggests that 
these observations are more consistent with a modulating region as large as 
35 a.u. with a larger diffusion coefficient than is usually assumed. 

2. Simple Model Including Time Dependent Diffusion . The basic features 

of the model are as follows 


(1) The modulation is treated in a one dimensional solar wind so that 
the basic coordinate is the distance X inside the boundary. 

(2) There are no energy loss effects (a consequence of one dimen- 
sional propagation) so that the behavior of the modulated differential density 
U(T) at a particular kinetic energy T can be analyzed independently of the 
behavior at other energies. 


(3) The modulation at a particular energy depends on three very 
sbwly varying parameters [i] the distance X from the boundary; [ii] the 
diffusion coefficient k which is taken to be independent of position inside 
the boundary for simplicity; [iii] the solar wind velocity V which is inde- 
pendent of position. 


Subject to these conditions the modulated density U(X,t) at time t 
obeys a time dependent Fokker-Planck equation of the form 


3U 3 2 U 9U 

at K ax 7 ax 


(i) 


Parker (1965) has analyzed the solution to equation (1) in terms of the 
probability density for finding a particle which crosses the boundary (X = 0 ) 
at t = 0 inside the boundary at X > 0 and t > 0. With a small change in 
notation we can write Parkers solution for a single particle in the diffu- 
sion approximation as: , j. v\2 




-(V t + x)' 
4<t 


( 2 ) 


where P is the probability density for finding a particle at X and t given 
that it is at the boundary at t = 0. Assuming a constant particle density 
U 0 at the boundary, the density observed at X and to is then 


U(X 


,to) = J to U 0 P(X,V,K,t 0 - t)dt 


(3) 
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where X, V, and < are all slowly varying functions of time. In the limit 
that the parameters are constant, the right hand side of equation (2) can 
be integrated with respect to t and equation (3) reduces to the familiar 
stationary solution U S (X) 


MX) = U Q e 


-vx 

K 


(4) 



The effect of time 
variations can be understood 
by examining the behavior of 
P (t) as illustrated in Fig. 

1 for two different values of 
k. P is plotted as a function 
of the time in the past at 
which a particle seen now (at 
t = t Q ) entered the modula- 
ting region. Note that the 
most probably time delay is 
about 34 days for k = 2 x 10 21 
cm 2 / sec and 40 days for 
< = 1.6 x 10 21 cm 2 /sec. 

Thus a 20% difference in k 
corresponds to about a 17% 
difference in time delay for 
the example shown. The im- 
portant point is that if 
the interplanetary medium 
characteristics (k,V,X) 
are changing in time, 
even slowly, particles of 
different < are sampling 
these characteristics at 

different times. Therefore the observed intensity over a range of energies 
(and k's) at a particular time corresponds to a set of different average 
modulating parameters occuring over a range of times in the past. This will 
result in a "phase lag" which depends on k(R) (and X and V). 

The effect is best represented formally in terms of the average time 
delay t which is somewhat longer than the most probable delay time indicated 
in Fig. 1. The average time is defined as follows: 


-100 -50 

* 0 ) (Doy«) 


Fig. 1. The distribution of delay times be- 
tween the detection of a particle at t = t D 
and its having entered the modulating region 
at some prior time t. The area under each 
curve is equal to e~VX/ic a nd the most prob- 
able time shifts to longer times for smaller 

K. 


T = 


f 


t P(t)dt 


f 


P(t)dt 


(5) 


It should be noted that changes in x produced by changes in P with time 
are second order effects and are to be neglected in evaluating equation 
(5). It is shown in Appendix I that to first order the solution of equa- 
tion (3) under slowly varying conditions at a given time t = t Q is given by 


U(X,t 0 ) = U s [X,t Q - t(k)] 


( 6 ) 
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At high energies k(R) becomes very large and x(T) becomes small and the ex- 
pected intensity is equivalent to the stationary solution. Equation 6 there- 
fore shows that time dependent diffusion provides a mechanism for producing 
a phase lag at low energies with respect to that at high energies, i.e. a 
"hysteresis effect". 

Equation (5) which 


(7) 

where the normalizing denominator in equation 5 has been explicitly written 
as e^/ K . Again for simplicity we consider only the limiting cases of large 
and small k. In the case of large < (i.e. VX/k << 4) the convection term in 
the exponent of the integral in equation 7 can be neglected in the integra- 
tion yielding 

n. y2 

t(X,k) (8) 

which is equivalent to the diffusive propagation time in a one dimensional 
array of fixed scattering centers. In the limit of small k (XV/ic >> 4) the 
diffusion term in the exponent of the integral in equation 7 is neglected 
in the integration with the result that 

t(X,V) =| (9) 

The limiting behavior predicted by equations (8) and (9) can be compared 
directly with observations to illustrate the concepts of the model. 

A detailed study of the rigidity dependence of the magnitude of hy- 
steresis effects for a fixed value of the Deep River neutron intensity has 
been compiled by Rygg, et al. (1973). These same data also have been analyzed 
to find the time lags between the proton and helium intensities and the neu- 
tron intensities which result in a single valued regression (i.e. a closed 
loop). The values derived from the observations are listed in Table I along 
with the corresponding rigidity, velocity and diffusion coefficient k (calcu- 
lated from the Mariner IV Interplanetary Magnetic Field Spectra as given by 
Jokipii and Coleman 1968). 


3. Evaluation of x and Comparison with Observations . 

can be written explicitly as: 

VX /•co / . 9 


t(X,V,k) = 


- (vt + X) 2 
4xt 


(4itk) 


TABLE 1 


Observed Hysteresis Phase Lags 



Energy 

Rigidity 

B 

x ca l(cm 2 /sec) 

x (days) 

Protons 

60MeV 

. 34 1GV 

.34 

1 x 10 21 

150 ± 30 


100-250 

. 600GV 

.54 

2.6 x 10 21 

150 ± 30 


260-720 

1.078GV 

.754 

4 x 10 21 

120 ± 45 


>720 

>1.368 

>.824 

>6.2 

45 ± 10 

Helium 

60MeV/nuc 

.682 

.34 

1.7 x 10 21 

150 ± 30 


100-150MeV /nuc 

1002 

.47 

2.8.x 10 21 

60 ± 15 


200-250MeV/nuc 

1376 

.59 

4.4 x 10 2L 

45 ± 15 


>260MeV/nuc 

>1492 

>.62 

>.53 

40 ± 10 
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These observations are 
compared in Fig. 2 with the pre- 
dictions of the limiting approxi- 
mations to equation 7 discussed 
above. The observed time lags 
are plotted as a function of the 
corresponding calculated value 
of k showing a gradual decrease 
with increasing k. The solid 
straight lines represent the 
diffusive progaga tion times 
(equation 8) in interplanetary 
regions of 10, 20 and 35 a.u. 
respectively. The curved dashed 
lines represent the limiting be- 
havior at small ic (equation 9) 
for these same regions. The ob- 
servations are seen to be crudely 
consistent with the areas bounded 
by the solid lines for 10 and 20 
a.u. However, the approximation 
for ignoring the convection term 
in equation 7 is breaking 

down in the region where the data 
are plotted. If the upper limit 
imposed by equation 9 is taken 
literally (it may in fact be wrong by factor ^2 because of our use of a 
one dimension model) , the largest observed phase lags of up to ^150 days 
imply a modulating region somewhat larger of 'v X = 35 a.u. Furthermore 
to be consistent with the model the characteristic diffusion coefficient for 
the observed particles over the whole region would be somewhat larger than 
the values plotted here. In particular if the average value of k for these 
particles were a factor of 'W larger, equation 8 would be quite valid and 
would lead to agreement with the observations in a region of 30-40 a.u. Such 
a result, modulation in large region with large diffusion coefficient, is 
consistent with current hysteresis observations and the small gradients 
currently being reported on Pioneer 10 (McKibben, et al. 1973). 

The model presented here provides an explanation of "hysteresis" 
which does not require that k be an inseparable function of R, t and V 
although such behavior may be part of the complete picture. Furthermore, 
more sophisticated analysis based on these concepts may provide an expla- 
nation for the anamolous behavior of the helium spectrum reported recently 
(Mason, Munoz and Simpson, 1973). 

Appendix 1 

To determine the effects of slowly varying interplanetary conditions 
we express the probability density P(t) as 

P(X,V,K,t) = P 0 (t) + 6P 0 (t) + ... A X ) 



•Fig. 2. 

Observed delay times for protons and 
helium nuclei during the 1969-1970 
Solar Maximum. The points are plotted 
at calculated values of k based on the work of 
Jokipii and Coleman (1968). The curves 
indicate the expected dependence of x 
on k for three positions of the modu- 
lating boundary under the approximation 
of the model. 


where 
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P Q (t) = P(X 0 , V , k d , t) is P 0 (t) evaluated under conditions existing at 
t = t Q and 


6 P 0 (t) 


3P n AX , SP^ AV_ , 6P^ Ax_ ( . 

3X At 0 + 3V At 0 3k At 0 { ~ 


A2) 


with AX/At Q , AV/At 0 and A</At 0 being the slow rates of change in each of 
these parameters. Under these conditions we can write equation 3 in the 
text as 

U(X,t 0 ) = jj U oP 0 (t 0 - t) dt U Q [t - t 0 ] Uo = 0 + ... A 3 ) 


The first term on the righthand side is simply U s (X,t Q ) by definition. 

We evaluate the second term subject to the condition that [l/P 0 * 6 P 0 /6 1 >>T 

as defined by equation (5). (In the stationary solution the assumption is 
that [ (1/P) ( 6 P /6 t Q ) ] -1 = Here for the first order contribution we simply 
assume it is very large but not infinite.) Thus we can take this derivative 
with respect to t 0 out of the integral and write the second term of A 3 as 



6 

C to 

u o 1 [t - t Q ] P 0 (t 0 - t)dt 


T U Q P 0 (t 0 - t)dt 

J-co 

J-CO 


a 4 ) 


where we have used the definition of t in equation (5) . Thus equation A^ 
can be rewritten as 

u(x,t Q ) = u s (x,t 0 ) - t ]^~ [u s (x,t 0 )I + ... = u s (x,t o - t) a 5 ) 

where we have evaluated <$/<5t 0 over a time interval x neglecting changes in 
T as second order effects. 
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Technical Report #73-132 
DIRECT OBSERVATIONS OF THE CHARGE STATES OF LOW ENERGY SOLAR PARTICLES 

G. Gloeckler 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland (USA) 20742 

C.Y. Fan 

Physics Department, University of Arizona, Tucson, Arizona (USA) 85721 

D. Hovestadt 

Max-Planck Institut fur Extraterrestrische Physik 
8046 Garching (Germany) 

The charge states of carbon and oxygen of solar origin have been 
measured directly in interplanetary space. At 100 keV per nucleon 
the C +5 /C"*" 6 and 0 +7 /0 +8 ratios are 1.8 and 1.6 respectively. We 
find the abundance ratios of low energy heavy nuclei to He signi- 
ficantly larger than corresponding photospheric values: The en- 

hancement of 0/He is 35 and both Si/He and Fe/He are overabundant 
by a factor of 100. To explain these observations we propose a 
mechanism which first preferentially accelerates heavy ions and 
is followed by either storage of these ions in the coronal regions 
or strong adiabatic decleration. 

1. Introduction . Chemical composition measurements of solar flare parti- 

cles above 10 MeV per nucleon, pioneered by the Goddard Emulsion Group (Fich- 
tel and Guss , 1961; Biswas and Fichtel, 1965), have been extended in recent 
years to very low energies using dielectric track recording devices (Fleischer 
and Hart 1973; Braddy _et al. 1973; Crawford j 2 t al. 1972). The most striking 
feature of these low energy (<10 MeV per nucleon) observations has been the 
enhancement of the ratios of heavy to light nuclei over corresponding photo- 
spheric ratios in contrast to results of the Goddard group at higher energies 
(>10 MeV per nucleon) , which showed that the solar particle composition re- 
flected photospheric abundances (Bertsch et al. 1973) . Mogro-Campero and 
Simpson (1972) find that 5-60 MeV per nucleon solar cosmic rays also show a 
heavy particle enrichment which increases with atomic number and varies from 
flare to flare. To explain the increasing enhancement of heavy nuclei both 
with decreasing energy and increasing charge various acceleration models 

have been proposed, all of which require an establishment of charge state 
equilibrium for ions in the acceleration process (Ginzburg and Syrovatskii 
1964, and references therein; Ramadurai 1971; Cartwright and Mogro-Campero 
1972). Direct measurements of the charge state of low energy ions are there- 
fore clearly important in order to check the validity of these models. 

In this paper we present our preliminary results on low energy charged 
particles observed during a small energetic particle event on October 17-18, 
1972. Our measurements provide the first direct evidence for partially 
stripped carbon and oxygen at 100 keV per nucleon and extend heavy ion mea- 
surements down to 50 keV per nucleon. 

2. Instrument Description. The data have been . obtained using the Univer- 

sity of Maryland Electrostatic Deflection Analyzer (Fan et_ al. 1971) on board 
the IMP 7 (Explorer 47) satellite which was launched September 22, 1972 into 
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a nearly circular, 240 000 km apogee orbit. Our method of particle identi- 
fication is based on the fact that the amount of deflection, d, of a non- 
relativistic ion with an effective charge Z and kinetic energy T in a known 
electrostatic field is given by d = gZ*/T. The constant g is determined by 
the geometry and voltage of the deflection system. 


In Fig. 1 we show the cross 
sectional sideview of the sensor 
designed to measure particles in 
the energy range 50 to 1000 keV 
per nucleon. The collimator de- 
fines the incoming direction of 
the ion, three separate electric 
field regions (two low field and 
one high field) provide the 
necessary deflection, and ten 
rectangular solid state detectors 
are used to measure the energy T 
and deflection d. A plastic 
scintillator cup is used to 
reduce penetrating charged par- 
ticle background. 

Fig. 2 illustrates the 
charge resolution of the instru- 
ment. The deflection vs. energy 
lines for each given charge Z* 
are broadened to reflect the 
energy and geometrical resolutions 
of the detector system. Dashed 
horizontal lines correspond to 
boundaries of sensitive areas 
of given solid state detectors 
and define the energy window for 
each detector. The electronics 
associated with the sensor iden- 
tify the triggered detector and 
record the pulse height of the 
signal proportional to the total 
energy of the particle at a maxi- 
mum rate of 0.9 events/sec. In 
addition, sectored counting rates 
corresponding to the absolute par- 
ticle fluxes in selected and fixed 
energy per charge windows are 
transmitted once every 5 sec. 

The pulse height informa- 
tion is essential for identifica- 
tion of the charge state of an 
ion. A projection of the shaded 
portions of deflection vs . energy 
curve in Fig. 2 for a given de- 
tector on the energy axis shows 


ELECTROSTATIC ENERGY - CHARGE ANALYZER 
UNIVERSITY OF MARYLANO 



Fig. 1. Cross sectional side view of 
the Electrostatic Energy-Charge Analyzer. 
The geometrical factor for each of the 
low field deflection regions is 3 mm 2 
sr. For the high field region it is 
1.8 mm 2 sr. 



Fig. 2. Deflection vs, energy curves for 
particles having charge states 1, 2, 6, 8 
and 26 in the high and low field deflection 
regions of the instrument. Proton and 
electron calibrations were used to derive 
these curves. Shaded areas show the allow- 
ed energy ranges in each detector and indi- 
cate the charge resolution of the sensor. 
The energy ranges are fixed by the position 
and geometry of the detectors , and the vol- 
tage of the deflection system. 
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that one can easily resolve charge states 1, 2, 5-8, and 20-26 over the en- 
tire energy range of each detector. In addition, by making reasonable assump- 
tions concerning relative abundances of heavy nuclei, it is possible to 
identify adjacent charge states in the 5-8 group in more restricted energy 
intervals . 


3. Results . Solar particles observed near earth on October 17 and 18 
were probably produced by a subflare on the eastern limb of the sun (S06, 
E66) beginning at Q830UT on October 16, 1972. A sudden commencement (1800 
UT October 18) and Forbush decrease (240QUT October 18) accompanied this 
event. The 10 minute averaged counting rate proportional to the flux of 
125-160 keV protons is shown in Fig. 3 and 
illustrates the complex structure of this 
event. The initial increase above back- 
ground occurs at about 0200UT on October 
17, some 17 hours after the onset of the 
parent flare, and is followed by two more 
prominent increases at 2000UT, October 17 
and 1100UT, October 18. In the insert we 
show that the satellite was outside the 
bow-shock region during the. major portion 
of the event. Hourly averages of the pro- 
ton to helium ratio computed at equal velo- 
cities (120 keV per nucleon) were found to 
remain constant at about 40 throughout the 
entire event. Time averaged abundances 
ratios should therefore be computed on an 
equal energy per nucleon basis. 



Our results on spectra and abundances 
of particles are based on data averaged 
over the 30 hour time period shown by the 
solid bar in Fig. 3. For reference we 
present in Fig. 4 the time averaged pro- 
ton and helium spectra. We note that these 


Fig. 3. Time profile for 125-160 keV 
protons for the three day period Oct- 
ober 17-19, 1972. The background rate 
of 0.65 counts/sec is most likely due 
to the sensitivity of solid state de- 
tectors to y rays generated in the 
spacecraft. The peak rate of 10 
counts/sec corresponds to a differ- 
ential flux of 1.3 x 10 4 protons/ 

(cm 2 sec sr MeV) , 


are best represented by a dependence on 

rigidity, P, of the form dJ/dP = J Q exp (~P/P 0 ) with different characteristic 
rigidities P Q . Since the proton spectrum is steeper than the helium -spectrum 
the p/a ratio depends on energy and goes from 14 at 20 MV to 2 at 50 MV. 


The spectra of heavy particles appear in Fig. 5. Assuming the nuclei 
to be fully stripped (Braddy, ert _al. 1973) we compute spectra for the CN0, 

Si and Fe group nuclei and show these, plotted as a function of kinetic energy 
per nucleon, along with the helium spectrum in Fig. 5(a). We find that the 
abundances of Fe and Si are about equal within errors of the measurement and 
that the Si/He and Fe/He ratios at 100 keV per nucleon are about 0.025, or 
larger by a factor of 100 than the corresponding solar photospheric values. 
Fig. 5(b) gives two point carbon and oxygen spectra computed under the assump- 
tion that these nuclei are fully stripped and that small amounts of B, N, 
and Ne are present. The C/0 value obtained is about 2 to 3 which is much 
higher than the usual value of 0.5 for solar cosmic rays (Teegarden et al. 
1973). Furthermore, using the dE/dx vs. E telescope ULET (Hovestadt and 
Volliaer 1971) flown on the same spacecraft we find for the identical time 
period the C/0 ratio at 500 keV per nucleon to be 0.5±0.1. 
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Fig. 5(c) shows the spectra 
of partially stripped nuclei. Re- 
taining the assumption that very 
little Be, B, N and Ne is present 
we have taken all Z = 5 to be CT* -5 
and all Z* = 8 to be 0 +8 . It is 
then possible to estimate the abun- 
dances of C+ 6 and 0 +7 . The (C** 5 + 
C f6 )/(0 +7 + 0 +8 ) ratio now becomes 
0.8 and is in reasonable agreement 
with measurements at 500 keV per 
nucleon. In Table I we list our 
results on abundances at 100 keV 
per nucleon along with previously 
published values at higher energies 
and photospheric abundances. A sys- 
tematic increase in the enhancement 
of heavy nuclei is evident both with 
decreasing energy and increasing 
atomic number of the particles. For 
example, 100 keV per nucleon solar 
flare 0, Si and Fe nuclei are over- 
abundant relative to the solar photo- 
sphere by factors of 35, 70 and 100 
respectively. 

4. Discussion. The enrichment 

of heavy low energy nuclei in solar 
cosmic rays relative to solar photo- 
spheric abundances , and the presence 
of significant amounts of partially 
stripped nuclei are indicative of the 
acceleration process and the physical 
conditions of the medium where the 
acceleration occurred. We can sug- 
gest the following model to provide 
a qualitative explanation for these 
observations . 



RIGIDITY P (MV) 


To produce an enhancement of . , , , .. 

r Fig. 4. Time averaged proton and helium 

heavy nuclex we require statistical differential rigidity spectra. To obtain 

acceleration in regions of sufficient these spectra we have corrected the pulse 

neutral hydrogen density for charge height distributions for background and 

equilibrium to be established (Ginz- used the fixed energy windows for each 

burg and Syrovatskii, 1964; Ramadurai detector to compute the flux. 

1971). Acceleration of this type can 
account for the observed energy de- 
pendence of the heavy particle enrichment. Under these conditions, the 
charge state of the accelerated particles will depend on their energy. For 
example, the equilibrium charge states of carbon and oxygen at 100 keV per 
nucleon are found to be 2.5 and 2.9 respectively. Our measurements on the 
other hand, indicate much higher values, specifically 5.4 for carbon and 
7.4 for oxygen. To account for this discrepancy we propose two alternatives. 
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(1) Accelerated ions traverse hot coronal regions where the proba- 
bility for electron pickup is vanishingly small, and are stripped of most of 
their remaining electrons. Given an integrated column density of about 
3 x 10 18 /cm 2 and a cross section for stripping of a K-electron of the order 
of 10- 2 1 cm 2 , we conclude that the ions must have been confined within the 
coronal region for about an hour before escape. 


(2) After acceleration and escape from the sun, the particles ■ are 
adiabatically decelerated (Jokipii, 1973). In this case, the measured charge 
states are the equilibrum charge states of particles at much higher energies. 
Our measurements of the abundance ratios C +8 /C+ 6 and 0 +7 /0 +8 would imply 
energy losses of the order of a factor of ten. 


(a) (b) (c) 



Table 1. Abundance Ratios for Solar Flare Particles at Different Energies 


Fig. 5. The time 
averaged differen- 
tial energy spec- 
tra for helium 
and heavy parti- 
cles. Data have 
been corrected 
for energy defects 
in the solid state 
detector and for 
background mea- 
sured independent- 
ly for each detec- 
tor with the de- 
flection voltages 
commanded off. The 
dashed and dotted 
curves represent 
a scaled helium 
spectrum and are 
drawn for refer- 
ence only, (a) 
Spectra of CNO, 

Si and Fe group 
nuclei assuming 
full stripping. 

(b) Carbon and 
Oxygen spectra. 

(c) Spectra of 
partially strip- 
ped nuclei. 


Abundance 

Ratio 


Solar Flare Particles 

Photosphere 8 

0.1 MeV/nuc. a 

1.5 MeV/nuc.* 3 

14-50 MeV/nuc. C 

C+5/C+6 

1.8 




0 +7 /0 +8 

1.6 




0/He 

0.25 

0.15 

0.025 

0.0068 

Si/He 

0.025 

0.012 

0.0025 

0.00035 

Fe/He 

0.025 

0.018 

0.0007 

0.00025 


^Present work. 

"Data compiled by Lanzerotti et_ a^. 1972. 
*“Teegarden et_ al . 1973. 

°Hange and Engvold 1970; Withbroe 1971. 
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The limited amount of data we have at present does not allow us to 
decide which of the two alternatives corresponds to actual conditions. 

More extensive observations of low energy ions will be required before de- 
cisive information about particle acceleration and their subsequent propa- 
gation to earth is obtained. 
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COHERENT PROPAGATION OF CHARGED PARTICLE 
BUNCHES IN RANDOM MAGNETIC FIELDS 

James A. Earl 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland (USA) 20742 

When the distribution function for energetic particles moving in 
interplanetary space is expressed in terms of eigenfunctions of 
an operator which describes pitch angle scattering by random 
magnetic fields, the familiar phenomenon of diffusion is the 
dominant solution of the transport equations provided that the 
spectrum of magnetic fluctuations is not too steep. However, 
if this condition is not satisfied, the dominant solution of 
the same equations is a new mode of transport in which density 
inhomogeneities propagate coherently along the guiding field 
with a characteristic velocity equal to half the particle speed. 
The most striking manifestation of this phenomenon is the ob- 
served "scatter free" propagation of solar flare electrons 
below 1 MeV. 


1. Introduction . For many years , it has been virtually taken for 

granted that cosmic ray propagation is fundamentally diffusive in nature. 
In early papers, the scattering centers responsible for diffusion were 
visualized as discrete magnetized clouds. However, in the physically 
more plausible configuration of random fields superimposed upon a smooth 
guiding field, diffusive transport, characterized by a finite diffusion 
coefficient, has also been established as a fundamental mode provided 
that small scale fluctuations among the random perturbing fields are 
sufficiently intense (Jokipii, 1966; Hasselman and Wibberenz, 1970; Earl, 
1973). More specifically, if the spectrum of magnetic energy density vs. 
wavenumber k which characterizes the perpendicular irregularities is a 
power law P(k) a k~^ with spectral index q, then diffusion occurs' when 
q < 2. This paper describes a novel mode of particle transport in which, 
for q ^ 2, density inhomogeneities propagate in a wavelike fashion with 
very little dispersion. Because this behavior is qualitatively differ- 
ent from diffusion, the new mode will be designated as "coherent". Since 
a complete treatment of these effects is available in a longer report 
with the same title as this one (U. of MD. , TR #73-098), only the most 
basic features of the coherent mode will be presented here with emphasis 
on physical interpretations rather than on mathematical proofs. 


2. Transport Equations . In the present one- dimensional case, the 

Boltzmann equation can be written as 


3f 3f 1 3_ < Ay 2 > 3f_ 

3t + yV 9z 2 3y At 3y 


( 1 ) 


where y is the cosine of the pitch angle relative to the guiding field, t 
is time, V is particle velocity, z is distance along the mean field, and 
f is particle density in phase space. (For momerttum p, this quantity is 
proportional to cosmic ray flux, J = f/p 2 .) The last term in Eq. 1 describes 
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pitch angle scattering as characterized by a Fokker-Planck coefficient 
given by Jokipii (1966) in the form, 

^ = A (1 - „*) Ml ' 1 (2) 

where q is the spectral index discussed earlier and A is a parameter which 
depends on the r.m.s. intensity of the random fields and upon particle mo- 
mentum and velocity. The key to understanding Eq. 1 is to recognize that 
the' scattering operator is such that Sturm Liouville theory is applicable. 
This ensures that there exist eigenfunctions R and relaxation times x 
which are defined by 


3_ 

3y 


< Ay 2 > 

At 


3R* 

3p 


R = 0 


(3) 


and that the distribution function can be expressed as an eigenfunction 


series f = 


Ef K R K 


in which the coefficients f K decrease in magnitude rapidly 


with < so that only the first few terms are needed to represent f quite 
accurately. Because the functions R K form an orthogonal set, it is easy 
to derive the following set of coupled differential equations for the co- 
efficients : 


Mr 

at 


= -v 01 


^ = - v 
at x 2 ] 


afi 

az 

< 

O 

CO 

3fq 

3z 

(4) 

lii 

at 

+ 

11 

-v 0 f 

Mq 

az 

-V 12 

af2 

az 

(5) 

af ! 

3z 

~ V 2 3 

3f 3 
3z 

(6) 

Ma 

at 

T 3 

-v 03 

.Mo. 

3z 

-V 2 3 

Mz 
a z 

(7) 


where Vqj, V 0 3 , V 12 and V 23 are characteristic velocities defined by inte- 
grals which would take too much space to give here but which are exactly 
analogous to quantum mechanical matrix elements. The coefficient f Q of 
the isotropic eigenfunction R 0 = const, is proportional to the spatial 
density U while the odd anisotropic components Rj , R 3 , ... give rise to 
a streaming flux S whose gradient is equal to the time rate of change of 
density. In diffusion, only the two lowest order components fo and f^ 
are important, but, in coherent transport, the components f 2 and f 3 also 
play a significant role. 


To make clear the difference between coherent 
and diffusive regimes, it is worth considering, for 
each case, how scattering causes the first order aniso- 
tropy Ri to change with time. In Fig. la, the Fokker- 
Planck coefficient defined for q < 2 by Eq. 2 (solid 
curve) dips to zero in the cusp at y = 0 , but the slope 
of Ri at this point is very steep. (See Fig. lb.) 
Consequently, a large derivative acts in the .scattering 
operator to compensate for weak scattering so as to 
produce the same fractional rate of change at y = 0 as 
everywhere else. Thus, according to Eq. 3, an aniso- 
tropy proportional to Rj will decay exponentially with 
no change in shape at a rate characterized by the re- 
laxation time Tj. As will be indicated later, this 
rapid decay of the first anisotropy is closely related 
to diffusion. In the coherent regime, q 1 2, the 
dashed curve describing (<Ay 2 >/At) not only has zero 
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value but also has zero slope at the origin. This means, in effect, that 
weak scattering occurs not just at a single point but over a finite region. 
Consequently, the compensation discussed above does not occur and the eigen- 
function Rj is the step function shown as a dashed curve in Fig. lb. More- 
over, the first relaxation time is infinite because the derivative of 
the step function is zero in regions where the Fokker-Planck coefficient 
is large, while near the origin where the derivative is large, the coeffi- 
cient vanishes. For q < 2, the situation is analogous to quantum mechan- 
ical tunneling in that particles are able to pass through the region of 
weak scattering at y = 0 to allow a first order anisotropy to relax. For 
q > 2 , this effect is inhibited so that an anisotropy proportional to the 
step function is stable. It is this persistence of anisotropy that gives 
rise to coherent effects. 


3. Diffusion . In standard works on transport theory (Weinberg and 

Wigner, 1958, p. 231), the diffusion approximation is obtained by neglecting 
the derivative (3f^/3t) as well as coefficients of higher order than fj. 
Thus, according to Eq. 5, the fj component reaches a secular equilibrium, 
fl = -T 1 V 0 i(3f o /8z) which can be substituted into Eq. 4 to give the diffu- 
sion equation 



Tl(Voi ) 2 




( 8 ) 


Evidently the diffusion coefficient D = TiCVqi ) 2 is well defined when 
is finite but it becomes infinite in the coherent regime. The solution 
of Eq. 8 for impulsive injection at t = 0 and z = 0 of a localized bunch 
of particles is the familiar Green's function for diffusion 

f Q = [F o /(4-rrDt) 2 ] exp {-z 2 /4Dt} (9) 

where F q is a measure of the total number of particles injected. 

4. The Coherent Regime. When q - 2 , the eigenfunctions R 2 and R 3 are 

not only important but also they are "degenerate" in the sense that their 
eigenvalues are equal so that x 2 = T 3 . Furthermore, V 03 = V 12 and = 00 . 
Under these circumstances, the symmetry of Eqs. 4 - 6 is such that the 
solutions must meet one of the following two conditions: 


f “ + f l J f 2 - + f 3 


or f 


= - ff; f 2 - - f; 


( 10 ) 


In the spirit of the diffusion approximation discussed above, the spatial 
and temporal derivatives of f 2 and f 3 can be neglected so that Eqs. 4 and 
6 combine to give an equation for f 


^ +V 01 III C = *3(V 0 3) 2 0 s2 - 

The transformation y = z - V 0 jt reduces this to 


(ID 



t 3 (yo 3 ) 2 


3 2 f 


3y 


+ 


9 2 f< ? + 


( 12 ) 
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which is of the same form as Eq. 8 . Consequently, the solution of Eqs. 11 
and 12 corresponding to impulsive injection can immediately be written by 
analogy with Eq. 9 as 

f Q + = [F o /(16TTD 5k t)^] exp {- (z - V 0 it) 2 /4D*t) } (13) 

where the diffusion coefficient has been replaced by a parameter D* = T 3 CV 03) 2 
that I call the coefficient of dispersion. In Eq. 13, the diffusive Gaus- 
sian of Eq. 9 centered at the stationary point z = 0 has been replaced by 
a moving coherent Gaussian centered at z - Vg^t. For isotropic injection, 
this localized pulse moving in the +z direction with velocity Vqi is accom- 
panied by a symmetrical pulse f~ moving in the -z direction. The latter 
disturbance satisfies the second condition of Eq. 10. 

The summation, implied by Eq. 10, of an isotropic component and a 
step function of equal amplitude gives for F*" a distribution function which 
is constant for p > 0 and zero for p < 0. (This form is slightly modified 
by f 2 and f 3 but these components are small compared to f Q and fj.) This 
describes a bunch of particles with velocities isotropically distributed 
over the forward hemisphere with no particles moving backward. The compo- 
nent f~ represents a similar group moving in the opposite direction with 
velocities in the other hemisphere. These bunches move coherently because 
the velocity vector of an individual particle wanders randomly over one 
hemisphere giving an average component along the z axis equal to that of 
other particles in the group. This speed is V 01 = (V/2). The steady in- 
crease of the width of the coherent disturbances, characterized in Eq. 13 
by the dispersion coefficient D*, results from statistical fluctuations in 
this averaging of directed random velocities in much the same way that or- 
dinary diffusion results from fluctuations in the averaging of nearly iso- 
tropic random velocities. 

To summarize, Fig. 2a shows schematically for the 
diffusive regime the profile of density vs distance z after 
the injection of a localized pulse at z = 0. Here, the 
narrow spikes shown moving outward with velocities ±V 01 on 
each side of the Gaussian arise when the time derivative 
(3fx/9t) is retained in a treatment similar to the simpli- 
fied calculation that led to Eq. 8 (Axford, 1965). The 
number of particles contained in these spikes decreases 
rapidly leaving the slowly evolving diffusive Gaussian as 
the dominant contribution. In the coherent regime depicted 
in Fig. 2b, the density in the coherent Gaussians decreases as 
they spread out, but the total number of particles in each 
is constant, and the diffusive Gaussian is completely ab- 
sent. The qualitative similarity between the coherent 
disturbances and the spikes of the diffusive regime suggests 
that, in the case where tj is finite but very large, a useful description 
can be obtained by simply substituting coherent Gaussians (Eq. 13) for the 
localized spikes. In this mixed regime, shown in Fig. 2c, the number of 
particles in the coherent pulses decreases as particles are scattered out 
of them into the diffusive profile. The mixed regime is of practical 


( a ) diffusive 



-V„t o + v 0 ,t 
Fif 2 
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importance in interplanetary space where direct measurements of an exponent 
q slightly less than 2 (Sari and Ness, 1969) imply such a large value of 
that coherent effects should be very significant. Moreover, recent 
theoretical work (Klimas and Sandri, 1972; Kaiser, Jones and Birmingham, 
1973) suggests that tj may be large but finite even for q > 2. In this 
case, coherent effects would be very prominent in the mixed regime, but 
purely coherent transport would not occur. 

5. Comparison with Solar Particle Observations. Within the present 

limitations on length, it is not possible to present and discuss the ex- 
plicit formulas which relate the parameters appearing above to the power 
spectrum of interplanetary magnetic fluctuations. However, an implication 
of these relationships is indicated in 
Fig. 3 which shows a calculated profile 
of the intensity of 38 KeV electrons at 
earth vs. time after injection at the 
sun. A magnetic power spectrum typical 
of those actually observed in space was 
assumed. In the mixed regime occuring 
here, the coherent pulse is followed 
by a slowly decaying tail which is, in 
essence, a diffusive effect. The fact 
that this profile corresponds in every 
detail to the phenomenon designated by 
Lin (1970) as "scatter free" propaga- 
tion of solar flare electrons makes it 
clear that the theory explains his ob- 
servations. However, his terminology 

seems inappropriate now because it is the existence of scattering within 
the forward hemisphere that ensures that all electrons have nearly the 
same average velocity parallel to the field and thus makes possible the 
propagation of coherent groups. 

Coherent propagation offers a con- 
vincing explanation for several other poorly 
understood aspects of solar particle behav- 
ior. These include, for example, the initial 
field-aligned anisotropy of solar particles 
recorded by neutron monitors (McCracken, 

1962) and the double peaked profiles re- 
corded in many events (Fichtel and McDonald, 

1967, Fig. 1). To demonstrate that coherent effects are not restricted to 
low rigidity electrons. Fig. 4 shows how a Gaussian given by Eq. 13 pro- 
vides a good fit to a flare profile observed by McCracken, Rao and Bukata 
(1967, Fig. 22) while their detector of 7.5 to 45 MeV protons was pointed 
toward the sun. The extreme anisotropy predicted above for the coherent 
pulse explains why the antisun intensity was relatively small. The diffu- 
sive tail observed after 0045UT, which is similar to that in Fig. 3, is 
also predicted. 

These interpretations demonstrate that coherent effects exist. How- 
ever, the real significance of this new understanding of solar phenomena 
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is to suggest that coherent propagation may also play an important role 
elsewhere in astrophysics. 
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IN SOLAR PARTICLE FLUX DURING JULY 1972 
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Department of Physics and Astronomy, University of Maryland 
College Park, Maryland (USA) 20742 

Solar energetic particles were observed with excellent time resolu- 
tion with the aid of a large hodoscope flown on a balloon. This 
paper reports on a transient increase in the intensity of solar pro- 
tons which can be interpreted in considerable detail as evidence for 
the coherent propagation of particles injected by a west limb flare 
occurring at 1025UT on July 22, 1972. 

1. Introduction. Balloons provide a combination of payload lifting capa- 

bility and exposure duration that opens up new possibilities for investiga- 
ting rapid temporal variations of primary cosmic ray intensity. This paper 
describes an experiment, designed to exploit these possibilities, which em- 
ploys a detector whose geometric factor is large enough to reveal changes on 
a scale of minutes in the quiet time flux of 100 to 400 MeV protons. On the 
one hand, payload weight limitations rule out comparable satellite instru- 
ments, while, on the other hand, this region of the primary spectrum is com- 
pletely inaccessible to large sea level detectors. Moreover, the balloon 
platform offers the advantage over satellites that intrinsic temporal varia- 
tions are not obscured by rapid changes resulting from orbital motion through 
the spatial structure associated with geomagnetic cutoffs. 

Recent theoretical work has stimulated interest in "cosmic-ray scin- 
tilla:ions" which are broad band intensity fluctuations arising, in much the 
same way as optical and radio scintillations , from random inhomogeneities in 
the geomagnetic and interplanetary magnetic fields (Owens and Jokipii, 1972) 
(Owens and Jokipii, 1973). The present experiment is well adapted to study 
these statistical effects, but this short report, which is based upon a pre- 
liminary analysis of the data, will focus, instead, upon an interesting short 
term effect associated with a flare injection of solar protons. Specifically, 
a brief increase (20 min duration) in the intensity of vl00 MeV protons was 
seen at a delay time consistent with that expected for the coherent propaga- 
tion of a localized group of particles from sun to earth (Earl, 1973). This 
interpretation suggests that the phenomenon designated by Lin (1970) as 
"scatter free propagation" occurs not only at the low rigidities of kilovolt 
electrons but also, on occasion, at high rigidities where only diffusive 
propagation has heretofore been considered. 

2. Instrumentation. The hodoscope is shown in Fig. 1. The detector ele- 
ments are 7 trays of 27 geiger tubes with successive layers perpendicular to 
one another, 3 plastic scintillation counters, and a lucite Cerenkov counter. 
Geiger tubes in the top four trays are individually connected to coincidence 
circuits which, when triggered by a particle whose trajectory is suitable 

for detailed analysis, initiate a readout cycle that takes 23.6 msec to trans- 
mit trajectory, pulse height and housekeeping information for each event. 
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Sandwiched between the 2nd and 3rd trays is 
a 2.5 cm thick plastic scintillator to mea- 
sure dE/dx for incident particles. Below the 
fourth tray, are two 14 cm thick plastic scin- 
tillators which record the ionization of 
stopping particles in two energy ranges; (90 
to 170 MeV/nucleon and 170 to 230 MeV/nucleon). 
Each is viewed by two RCA 8055 photomultiplier 
tubes. Between these calorimeter counters, 
is a 5th tray (not shown in Fig. 1) to define 
those events which stop in the first counter. 
Two additional crossed trays below the second 
scintillator define the upper limit of 230 
MeV/nucleon for the detection of stopping 
particles. Analysis, however, proceeds in 
two additional energy ranges defined by the 



presence or absence of an output from a 5 cm 

thick Lucite Cerenkov counter viewed by nine The hodoscope detector, 

photomultiplier tubes. The threshold of this detector is 400 MeV/nucleon. 


Many details of this system are similar to those of smaller instruments 
flown in the past which incorporated 15 cm square detector elements with 9 
tubes in each tray (Rygg and Earl, 1971; Earl, Neely and Rygg, 1972). How- 
ever, the geometrical factor of 200 cm 2 sterad resulting from the present in- 
creased size of elements (50 cm) is two orders of magnitude larger than typi- 
cal satellite detectors and 10 to 20 times that of balloon-borne detectors 
traditionally employed to detect protons and helium in the energy range of 
interest here. The telemetry system, operated at a 16 KHz bit rate, provides 
nearly error-free transmission of approximately 30 events/sec. with multi- 
parameter pulse height analysis and trajectory information available for each 
event . 


Representative flight data from a one 
hour period at ceiling are presented in Fig. 
2 as contour lines plotted on a plane whose 
ordinate and asscissa are pulse heights in 
the dE/dx and calorimeter #2 counters, re- 
spectively. The numbers on these axes re- 
fer to logarithmic intervals each of which 
represents a factor of (4/3). Thus, 8 
channels corresponds to a factor of 10. 

The number of events per cell associated 
with each curve refers to particles which 
stopped in calorimeter #2 with no evidence 
of multiplicity in any of the trays. The 
clear resolution of components evident in 
Fig. 2 is, in part, due to the discrimina- 
tion against interactions made possible 
by the hodoscope, and, in part, due to the 



02 46 8K) 12 14 16 

CALORIMETER NUMBER 2 


Fig. 2. Pulse height matrix for 
particles stopping in the second 
calorimeter counter. 


good resolution in ionization measurements made possible by thick scintillators 
and efficient light collection. The intensities presented later of 90-108 MeV 
protons were derived from matrices, similar to that in Fig. 2, in which were 
collected particles stopping in the first calorimeter. 
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3_j The Coherent Event. The instrument has been flown on balloons 

launched at Churchill, Manitoba on the following four dates: July 18, 1971; 
July 21, 1971; July 21, 1972; July 26, 1972. All of these flights were suc- 
cessful. Both 1972 flights recorded solar particles from the active region 
that later gave rise to the great events of August 1972, but the phenomenon 
to be discussed here was probably not associated with that region. 


The coherent event occured 
during Flight 1330 on July 22, 1972. 
Starting at about 0500UT , a slowly 
varying influx of solar particles be- 
gan. This influx, which probably 
came from a flare on the back of the 
sun, shows velocity dispersion and 
is worthy of study on its own merits, 
but the feature of real interest ap- 
pears in Fig. 3 where the flux of 90- 
108 MeV protons increased briefly by 
about a factor of 2 over the nearly 
constant flux of solar protons from 
the earlier event. Both the onset 
of this increase at 1150UT and the 
termination at 1109UT took place in 
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Pig. 3. Solar proton fluxes on 
July 22, 1972. 


less than one minute which is the interval during which counts were accumula- 
ted for the points in Fig. 3. Because these abrupt changes seemed remarkable, 
we looked into the interplanetary configuration at the time of the event and 
found the situation depicted in Fig. 4. There was a small flare at 1025UT 
almost exactly at the foot of the interplanetary 
magnetic field line which passed through earth. 

This flare emitted x rays detected by an ionization 
chamber flown as an independent experiment on the 
same balloon (H.R. Anderson, private communication). 

By coincidence, the asymptotic viewing direction 
from Churchill for 100 MeV protons (McCracken, 1962, 

Fig. 1.2) happened to be within 10° of the field 
direction. 
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Fig. 4. The solar neighbor- 
hood during the coherent 
event. 


The facts that the increase in Fig. 3 was 
most prominent in the lowest energy channel as is 
expected for a steep solar particle spectrum and 
that it occurred shortly after a well documented 
flare leave little reason to doubt the association 
between the event and the west limb flare. However, 
the temporal profile of the pulse is completely 
incompatible with diffusion models that are frequently invoked in the interpre- 
tation of solar particle injections. On the other hand, the qualitative form 
and quantitative details of the profile are predicted as manifestations of a 
novel coherent mode of propagation (Earl, 1973). In brief, particles can propa- 
gate from sun to earth along interplanetary magnetic field lines in localized 
groups which move with characteristic speeds comparable to the particle velocity. 
In particular, it is predicted that the leading edge of the coherent disturbance 
will move with speed 0.8V = 0.34C where V. = 0.42C is the velocity of a 100 MeV 
proton and C is the speed of light. This speed is in good agreement with a value 
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of 0.36C calculated from the 1500 sec interval between flare and pulse onset 
assuming a typical distance of 1.1 AU between sun and earth along the inter- 
planetary field line. The pulse in Fig. 3 does not resemble the characteris- 
tic Gaussian profile expected for a coherent disturbance long after injection. 
However, the observed square profile is consistent with the shape predicted 
shortly after injection under conditions such that the intensity of interplane- 
tary magnetic fluctuations is very low. The angular distribution of particles 
in the coherent disturbance is very anisotropic with maximum flux being ex- 
pected along the field line direction coming from the sun and zero intensity 
being expected from the antisun direction. Fortunately, our viewing direction 
(See Fig. 4) was just right for observing the coherent effect but we have no 
way of telling whether or not the intensity was small in the opposite direction. 

4_. Discussion. The effect noted here is qualitatively very similar to the 

so called "scatter free propagation" of low energy solar flare electrons (Lin, 
.1970). Moreover, the theory which describes coherent propagation not only gives 
a convincing explanation of both these phenomena but also it relates them to 
conditions in interplanetary space. In particular, the spectral index q £ 2 
of the magnetic power spectra observed by Sari and Ness (1969) is such that 
prominent coherent effects are expected, and the pronounced variations observed 
in magnetic fluctuations (Siscoe, Davis, Coleman, Smith and Jones, 1968) sug- 
gest that, on some occasions such as the present one, the spectral intensity 
is so low that nearly pure coherent propagation could occur. However, we have 
access to no magnetic data that would confirm this suggestion. In any case, 
the identification of a coherent effect at relatively high rigidities should 
stimulate further work leading to improved understanding of the acceleration 
of energetic particles in solar flares. In this connection, it is worth noting 
in Fig. 4 that the Pioneer 10 spacecraft was quite close to the same line of 
force on which the coherent pulse was seen at earth. Thus, it might prove 
fruitful to search for a similar pulse recorded there after a delay determined 
by particle velocities and the distance between spacecraft and earth. 

Solar particle influxes which arrive at earth via the coherent mode are 
much more closely related to phenomena on the sun than are those which arrive 
after diffusion. Subject to the validity of the interpretations given above, 
we exploit this fact to conclude from the abrupt onset in Fig. 3 that the in- 
jection, at the sun, of 100 MeV protons occurs on a time scale of less than 
one minute. 
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ABSTRACT 


In an attempt to reduce divergences among independent measurements of the 
electron spectrum above 10 GeV, electrons were detected with the aid of a 
large hodoscope which incorporates 8 trays of Geiger tubes and 10 scintilla- 
tion counters into a calorimeter containing 33 radiation lengths of lead. 

This instrument has been thoroughly calibrated up to 9 GeV at the Cornell 
Electron Synchrotron and has been flown twice at a depth of 5 g/cm 2 for a 
total of 30 hours at ceiling. Electron showers at 10 GeV are completely 
contained within the calorimeter where their ionization is sampled twelve 
times, and where hodoscope trays provide information on the geometrical 
distribution of shower particles. Reliable identification of electrons was 
ensured by extensive information on individual events together with good 
statistical accuracy resulting from the long exposure time and the large 
geometric factor (335 cm 2 sterad). Electrons appeared in a histogram of 
number of events vs. ax 2 measure of goodness of fit to electron shower 
curves as a clearly resolved peak standing well above the nuclear background. 
Detailed characteristics of events in this peak are very similar to those 
of electron events recorded during calibrations. In particular, the ini- 
tial ionization and the lateral spread of shower particles are consistent 
with those of calibration electrons. In contrast, the behavior of the nu- 
clear continuum associated with minimum ionizing incident particles is dif- 
ferent from that of calibration electrons but similar to that of nuclear 
events associated with incident helium nuclei. The differential electron 
spectrum from 10 to 100 GeV was obtained and will be presented. 
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CHEMICAL ABUNDANCES OF COSMIC RAYS >4.5 GV MEASURED WITH A LARGE 
AREA PROPORTIONAL COUNTER-SCINTILLATION COUNTER STACK 

Jacques L’Heureux and C.Y. Fan 

Department of Physics, University of Arizona 
Tucson, Arizona (USA) 85721 

G. Gloeckler* and R. Mainardi 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland (USA) 20742 

A 6500 cm^-ster cosmic ray detector consisting of twelve gas counter 
trays sandwiched between two large area circular scintillation 
counters was flown from Palestine, Texas in November of 1972 to study 
the composition of primary particles >1.5 GeV/nucleon in the charge 
range 3 to 30. For each analyzed event, we record i) the particle 
trajectory using four 20 wire proportional counter trays, ii) the 
energy loss in each of the solid counters and iii) the dE/dx losses 
in each of the twelve gas counters. The large dynamic range of the 
detector is established by operating six of the gas counters in the 
ionization mode. A description of the instrument and some preliminary 
results are given. 

1^ Introduction . Present cosmic ray propagation models (see e.g. Shapiro 

and Silberberg, 1970) suggest that the measured abundances of C,0, and Fe at 
the earth are closely related to their compositions at the sources , whereas 
the amount of lighter elements (Li, Be, and B) , which are the fragmentation 
products of the heavier ones, reflects how much material was traversed by the 
primaries before they reach the earth. The composition versus energy for 
particles below 10 GeV/nucleon has been quite well studied and the results 
show a similar spectral shape for all elements. For particles of energies 
above 10 GeV/nucleon, however, the picture is quite different. 

Many groups have in the past year reported steepening differential 
energy spectra for secondary as compared to nuclei of mostly primary origin 
like C,0, and Fe (Balasubrahmanyan and Ormes , 1972; Juliusson et al. , 1972; 
Webber et al. , 1972; Smith et al. , 1973). Some of these groups have in addi- 
tion reported a slightly flatter Fe spectrum when compared to C and 0 (Bal- 
asubrahmanyan and Ormes, 1972; Webber et al. , 1972; Juliusson and Meyer, 1973) 
It is therefore clear that important information regarding particle sources 
as well as particle propagation in the interstellar medium can be obtained 
from high energy composition measurements . 

The instrument which we have designed to study this problem has the 
following desirable characteristics: 1) a large geometrical factor to detect 

a sufficient number of relativistic heavy cosmic, rays and ii) a low mass 
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to minimize the fragmentation probability of the incident particles in the 
instrument and to reduce the unwanted background induced by particles which 
may enter from the side. The light weight is also necessary to permit the 
instrument to reach high balloon altitudes. 


The instrument was constructed and flown on a high altitude balloon 
on November 4, 1972. The analysis of the data is in progress. This paper 
reports the preliminary results. 


2 . 


Ins trumentation. 


A schematic diagram of 
the instrument is shown 
in fig. 1. The main tele- 
scope consists of two cir- 
cular plastic scintilla- 
tion counters (NE-102) 
labelled as SI and S2. 

Each one is set in a 
light diffusion box and 
is viewed by four photo- 
multiplier tubes. Each 
scintillator is used for 
a dE/dx measurement of 
the incoming particles. 

The geometry factor is 
6480 cm 2 -ster. Sand- 
wiched between these 
counters , there is a 
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Fig. 1. Schematic cross section of the instrument. 


stack of gas proportional/ionization counters. The gas counters are in the 
form of 12 trays, each tray being subdivided into 20 separate counters. The 
square trays are 102 cm on the side and 5.1 cm thick. The 20 individual cham- 
bers are identical in dimension, 5.1 by 5.1 by 102 cm, each of which has a 
single high voltage wire at the center. The ground planes between trays and 
the partitions between the counters are made of wires spaced 1.27 cm from 
each other. All of the trays fit inside a common gas enclosure filled with 
a mixture of 90% argon and 10% methane at a pressure of one atmosphere. The 
gas pressure is maintained to an accuracy of 0.25%. 


The twelve trays, although identical in construction, are assembled 
with their center and partition wires alternatively in the X and Y directions. 
They are also operated in one of two modes. The top two and the bottom four 
are operated at 1550 volts in the proportional mode while the remaining six 
are in the ionization mode at 550 volts. The ionization trays are used for 
the dE/dx measurement of VH and WH cosmic ray nuclei (Z>20) since they are 
not expected to saturate at high ionization levels. Four of the six propor- 
tional trays (the top and the bottom two) are used to locate the X-Y position 
of the particle trajectory. This position information is used to compensate 
for path length differences, to correct for any non-uniformities in the scin- 
tillation and in the gas counters. The outputs from each individual wire 
chamber are summed for each tray and each of the 12 resultant pulses is digi- 
tized with a 10 bit A-D converter. 
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For each particle which penetrates through the detector system, we 
perform two dE/dx measurements in the plastic scintillation counters and six 
dE/dx measurements in the gas counters. For particles of charge less than 15, 
the six proportional counters are used while for nuclei of higher charge, an 
additional six measurements are made in the ionization chambers. Since the 
dE/dx in a solid scintillator is relatively constant over the entire relativistic 
range, the two dE/dx measurements are sufficient for the charge determination of 
the particle. For the energy determination, we attempt to use the known fact 
that the ionization loss in a gas counter by a relativistic particle increases 
logarithmically with energy in the range from 2 to 200 GeV/nucleon (Crispin 
and Fowler, 1970) to measure their spectra in this limited range. 

The first flight with this instrument took place on November 4, 1972 
from Palestine, Texas. The experiment was aloft for a period of 7.5 hours, 

5.7 of which were at an altitude of about 129,000 ft (3.2 g/cm 2 ) . The verti- 
cal cutoff rigidity at the gondola position varied from 4.5 to 3.1 GV. 

3. Data analysis and results. 

During the 5.7 hours at the floating 
altitude, 120,000 events of charge 
greater than 3 were analyzed, 

32,000 of which had a clean tra- 
jectory and were used for the deter- 
mination of the composition spectrum. 

The non-uniformities in the scintil- 
lator response vs. position were 
studied, using in-flight Carbon and 
Oxygen nuclei. The response function 
thus obtained was used in subsequent 
corrections to the energy loss 
measurements. A charge spectrum, 
obtained using only the scintilla- 
tor outputs corrected for the non- 
uniformities and for the path 
length differences , is shown in 
Fig. 2. In this plot. Carbon has 
a full width at half maximum 
(FWHM) of 7.6% corresponding to 
±0.2 charge units. The resolution, 
as expected, varies, as 1/Z as we go 

to nuclei having a higher charge. The saturation effect in scintillation 
counters however degrades the resolution reducing it to ±0.5 charge units at 
Fe. Our preliminary charge composition is in agreement with earlier work. 

The saturation effect in NE-102 scintillation material is shown in 
Fig. 3, where we have plotted the ratio of the light output to the energy 
loss (dL/dE) , as a function of energy loss. The energy loss scale is relative 
to that produced by a singly charged relativistic particle. Notice that the 
light reduction at Fe is 0.42, crowding the elements by a factor of 2.5 in 
this charge region. 
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Fig. 2. Charge spectrum from the scin- 
tillators corrected for their non-uniform 
response. 
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The analysis of the wire 
chambers for the energy measure- 
ment is yet to be completed. To 
illustrate the feasibility, a 
scatter plot of the scintillator 
dE/dx vs. the gas counter dE/dx 
is shown in Fig. 4. The gas 
counter dE/dx is the average of 
two wire trays which seemed to 
give the best performance during 
the flight. Clusters corres- 
ponding to the elements from Be 
to Ne are marked. The solid 
curves in the plot indicate the 
expected responses vs. energy 
when the logarithmic relativis- 
tic rise in the energy loss in 
the gas chamber is included. 

We hope that the results will 
show a marked improvement when 
all six gas trays are used. 
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Solar flare particles in the energy range from 125 to 160 KeV 
per charge were detected and analyzed by an electrostatic 
deflection vs. energy instrument on board the Explorer 47 (IMP 7) 
satellite. A harmonic analysis of the counting rates shows that 
frequently, the fluctuation of the particle flux can be 
approximately represented by N(t) = f(t)sin 2rrvt where f(t) is a 
smoothly varying function of t. The oscillation frequency v 
varies from flare to flare and is of the order of 2 x 10 -J ~/sec. 
The physical implications of the finding will be discussed. 

JL Introduction . The propagation of solar flare particles has been 

studied both theoretically and experimentally for many years. It is now 
generally agreed that, after the impulsive injection of the particles from 
the sun, the motion of the particles in the interplanetary space can be de- 
scribed as the superposition of a diffusion in a static magnetic field and 
a convection with the field lines carried outwards by the solar wind. The 
diffusion tensor K can then be related to the power spectrum of the inter- 
planetary magnetic field, (see e.g. J. R. Jokipii, 1971). In these studies, 
however, no attempt was made to investigate transient phenomena such as 
intensity variation due to solar wind fluctuations, shock waves, and hydro- 
magnetic waves, although these phenomena may contain important information 
relating to the mechanism of releasing particles from the active regions of 
the sun, or the modulation of particle intensity by magnetohydrodynamical 
waves. The purpose of this paper is to report a new phenomenon of the in- 
tensity variation of solar flare particles observed with a detector system 
on the IMP 7 satellite. It was found that frequently the intensity showed 
sinusoidal fluctuations of frequencies x 10"^~/sec, and the durations 
of the fluctuations were about one hour. Ihe physical implication of the 
finding is discussed. 

2 . The Detector System and the Satellite . The detector system used 

for the experiment is an electrostatic deflection vs. energy instrument. 

A schematic drawing of the instrument is shown in Fig. 1. It was designed 
to detect ions, electrons and positrons with energies as low as 100 KeV/ 
charge. Particles which enter the collimator are deflected in one of the 
three electrostatic field regions and strike one of the ten detectors. The 
amount of deflection d of a particle in the field regions is related to the 
energy T and the charge Q of the particle by the following expression: 
d = gQV/T, where V is the voltage difference between two deflection plates 


30 



2 

416 



of the given field region, and g is the characteristic constant. Therefore, 
iri this instrument, the range of the energy per charge of an ion, A(T/Q), 
is determined by the position and the size of the detector triggered. T 
and Q can thus be determined separately by an additional measurement of the 
energy deposit of the ion in the detector. The details of the system have 
been published elsewhere (Fan, Gloeckler and Turns, 1971). 

The voltages applied to the deflection plates are -10 KV, +10 KV. 

At these voltages, the ranges of the energy per charge of positive ions and 
electrons for the ten detector channels are from 0.125 to 1.90 MeV/Q. How- 
ever, for the solar flare events reported in this paper, only the 0.125-0.160 
MeV/Q energy channel, Pi, recorded rates high enough for the study. 

The apparatus was launched into space on September 22, 1972, on 
the IMP 7 satellite whose orbit has an apogee 38 earth radii and a perigee 
34 earth radii. Because of its large near circular orbit, solar flare par- 
ticles can be measured even when the satellite is in the tail region of the 
earth's magnetosphere. 

3. The Data, Method of Analysis and Results . The counts of the Pi 

are accumulated over four readout periods, i.e., 4 x 1.27843 sec, for the 
transmission by the satellite data system to the ground. These counts, 
after the subtraction of quiet time counting levels, yield the intensity 
profiles of solar events in these channels. Fig. 2 shows the Pi rate for a 
portion of the event observed on day 287 (October 13), and Fig. 3 gives the 
Pi rate observed on day 292 (October 17). The positions of the satellite 
in these time intervals are indicated by the sun-earth-probe angles, desig- 
nated as [_ SEP in the Figures. The rates in higher energy channels are 
much lower and are not presented here. The rates in the electron channels 
are practically the same as the background level. From the pulse height 
information, it was found that more than 95% of the particles recorded are 
protons. 


One distinguishing feature of the intensity- time profiles of 
these events is that they appear to be very rugged, having an appearance 
of saw teeth. The amplitudes of the fluctuations are far larger than 
statistical. In order to understand the nature of these features, an auto- 
correction function is constructed for the counting rates in each of the 
peaks shown in Fig. 2 and Fig. 3. Let N(k”r) = counting rate recorded at 
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t = kr, where k is a positive integer and t is the time interval between 
two neighboring data points. The quantity t is usually an integer m times 
the basic readout period t = 1.27843 sec. If K = total number of data 
points, the correlation function C is defined by the following expression: 


K-Ak-1 

C(Ak) = £ N(kr)N((k + Ak)r). 

. k=0 

TWo selected samples of the function are displayed in Fig. 4. These results 



Fig. 4. Correlation functions of the counting rates plotted in Fig. 2 
and Fig. 3, respectively. 


indicate that in these two time intervals, the Pi counting rates can be 
represented by N(t) = f(t)sin 2rrvt where f(t) is a smoothly varying function 
of t, and v's for these two events are respectively 2.45 x 1CT 3 and 1.66 x 
10~ 3 ~/sec. We have examined six other similar cases which were observed 
in the period from September 30 to November 6, 1972 (Days 274-311) for 
which the data tapes were available to us. The frequencies of the oscilla- 
tion varied from 1.5 x 10 -3 - 3.0 x 10 3 ~/sec. The physical implications 
of the oscillation will be discussed in the next section. 

4. Discussion . The intensity oscillations which we have observed 

cannot be due to the release mechanism of the particles at the sun. This 
conclusion is based on the fact that the energy window of the Pi channel is 
<v35 keV, and thus the dispersion in the arrival time for the particles in 
this energy band from the sun to the earth's orbit will completely smear 
out any possible oscillations whose periods are only of 5-10 minutes. We 
have also examined the nature of modulation of an anisotropic flux of 
particles by the spin of the spacecraft. We averaged the counting rate 
over m spacecraft revolutions, and were able to show that the oscillation 
amplitude of the correlation function due to the spin modulation for a flux 
of a degree of anisotropy "A" is ~A 2 /(mTT ) 2 with a frequency v b equal to the 
beat of the readout frequency (= 1/1.27843) "and spin frequency of the space- 
craft (1/1.3066). Because the beat frequency is 1.6 x lO’^'sec, about 
five times higher than the observed ones, and the observed amplitude is 
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Independent of m and is one order of magnitude greater than A2/(mrr)^ for 
m *= 16 and A = 1, we discarded this explanation. Therefore, it is more 
likely that the observed oscillation is due to the modulation of the solar 
flare particles by the magnetic fields in the vicinity of the spacecraft. 


Large amplitude fluctuation in the magnetic fields in the earth's 
magnetosheath was first reported by Siscoe et al. (1967). It was suggested 
that magnetic oscillations in the transition region are excited throughout 
the hydromagnetic frequency regime by the interaction between the solar 
wind and the magnetosphere. At the time of the measurement, the frequencies 
near 5 x 10 -3 ~/sec showed the greatest activity. From the fact that (1) 
all the particle intensity oscillations reported in this paper were detected 
during the time when the spacecraft was in the tail region of the magneto- 
sphere or near the bow shock and (2) the frequencies were nearly equal to 
5' x 10 --'/sec, we conclude that the hydromagnetic waves and the particle 

intensity oscillations are two closely related phenomena (Jokipii, 1973). 

The actual relation between the two phenomena is unknown at pre- 
sent. However, we can suggest the following two possibilities: 


(a) The energy spectrum of the particles was modulated by the 
waves. If the spectrum was as steep as what was observed on day 287, 
(shown in Fig. 5) , a small degree of 


acceleration or deceleration in the 
particle energies by the hydromag- 
netic waves would manifest itself as 
a large amplitude oscillation in the 
particle flux recorded by the 
detectors. 

(b) The particles were 
confined to move along bundles of 
magnetic lines of force but there 
was a non-uniformity in the particle 
density. An intensity oscillation 
would then be a natural result of a 
"flapping" motion of the lines of 
force due to the hydromagnetic 
waves. 


Fig. 5. Energy spectrum of 
solar flare proton. 
He, and CNO ob- 
served on Day 287, 
1972. 



Mev/nucleon 
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Heavy interplanetary particles (3<Z<30) have been 
observed during a solar flare period between October 29 
and November *1, 1972. The nuclear abundance has been 
determined by a two parameter analysis (dE/dX vs E 
technique) of a thin window proportional counter - solid 
state detector device aboard the Explorer 47 (IMP-7) 
satellite. Energy spectra are given for C, 0 and Fe. 

The elemental abundance relative to oxygen is deter- 
mined for carbon, neon, magnesium, silicon and the 
iron group in the energy range 0.62 to 6.9 MeV/nucleon. 

For the Fe/0 ratio no energy dependence could be observed 
within the error limits. 

1. Introduction. For several particle flares attempts have 

been made in the past to determine energy spectra and relative 
abundances of solar cosmic rays. Emulsion techniques at high 
energies (above 20 MeV/nucleon) ( Bertsch et al. . 1973a and 
references therein), active satellite experiments at energies 
down to about 8 MeV/nucleon ( Mogro-Compero et al. , 1972a, Tee- 
garden et al. , 1973), and etching techniques at energies down 
to 100 "keV/riucleon ( Price et al. , 1971, Crawford et al. , 1972, 
Fleischer et al. , 1973) revealed a rough agreement of the cosmic 
ray abundances with the solar abundances. A considerable over- 
abundance of iron in solar cosmic rays has been observed at low 
energies. Crawford et al. . 1972 and Braddy et al. , 1973 report 
an energy dependent increase of the Fe/Si and Fe/0 ratio in the 
energy region 0.30 to 20 MeV/nucleon, Measurements of Mogro- 
Compero et al. , 1972b indicate a high variability in the Fe/0 
ratio from flare to flare above 5 MeV/nucleon. 

The present paper for the first time reports abundance 
measurements obtained with an active satellite experiment which 
coversthe low energy range down to about 90 keV/nucleon (for 
iron) The results presented in this paper go down in energy only 
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to about 600 keV/nucleon. Results on lower energies will be re- 
ported later, since the use of the dE/dx data of Northcliff et 
al. , 1970 below this energy becomes increasingly questionable 
and a detained analysis of calibration measurements will be 
necessary. 

2. • Instrumentation and satellite. The experiment is a dE/dx 
residual energy telescope using the proportional counter as a 
dE/dx device and the surface barrier detector for determination 
of the residual energy. A plastic scintillator cup serves to 
reject penetrating particles. The lowest energy limits for each 
nuclear species is determined by the total thickness of material 
in front of the solid state detector, which amounts to about 
0.328 ug/cm^ polyethylene equivalent, and by the electronic 

2 

thresholds. The effective geometrical factor is 1.03 cm ster. 

The isobutane gas density of the proportional counter was 
stabilized to within-, less than 0.5% over the two months of 
operation after launch .of IMP-7. For more details see Hovestadt 

et al. . 1971 . 

From calibration measurements and response calculations it 
turns out that at low energies, where the heavy particles are 
not fully stripped any more, the response tracks in a two 
dimensional dE/dx - E matrix merge or cross each other for 
different nuclear charges. Nevertheless at entrance energies 
above about 400 keV/nucleon a clear separation of individual 
even Z nuclei is possible up to iron. For higher energies the 
charge resolution is mainly determined by path length variation 
due to the finite opening angle of 21 degrees half-angle. 

The satellite was launched on September 21, 1972 into a 
nearly circular orbit at a radial distance of about 38 earth 
radii. Because of the orientation of the instrument in the space- 
craft and the satellite attitude the telescope scanned all 
directions within the ecliptic. ' 

5. Observational result. The data were obtained during a 
weak solar particle flare in the period October 29 to November 4, 
1972. A total of 24470 events have been accumulated over the 
entire period of increased heavy particle flux as indicated on 
figure 1 by a bar. The flux increase was very likely associated 
with an optical flare region near the central meridian of the 
sun which occurred on October 29, 1972 at about 16:00 UT. The 
first increase of heavy particle flux was observed about 9 hours 
after the optical flare. Figure 2 shows the two dimensional 
pulse height distribution with D1 representing the pulse height 
in the proportional counter and D2 in the silicon detector. 

Tracks for carbon, oxygen, neon, silicon, magnesium and the 
iron group are clearly separated. 
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0 t2 0 12 . 0 12 0 12 0 12 0 12 0 12 Tim« 0 

303 (Oct. 29 ) 304 305 306 307 308 309 


Fig. 1. Heavy ion counting rate (Dl-single rate) versus time 
(response: carbon 0.15 - 4.6 MeV/nucleon; 0: 0.12-9 MeV/nucleon; 
Fe:>0.09 MeV/nucl.). The position of the satellite with respect 
to the sun-earth line is indicated by SEV. 


# r-ir* 



Fig. 2. Two dimensional dE/dx versus E pulse height distribution 
of the flare October 29 to November 4, 1972. Total number of 
event = 24417. D1 = proportional counter, D2 - silicon detector. 
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The number of particles along each track has been converted 
into ENERGY SPECTRA by using a calibration curve of incident 
energy versus D2 pulse height, based on Northcliff’s et al. , 

1970 dE/dx-tables for heavy ions, and taking into account pulse 
height defects of the silicon detector as determined experi- 
mentally at a heavy ion accelerator. 


Figure 3 shows the 
differential energy spectra 
for carbon, oxygen and iron. 
Along with the spectra is 
drawn a curve with the form 

E”^, which corresponds to the 
spectral form of all three 
spectra. 

Within a wide energy 
band as indicated in table 1 
for each species the 
ABUNDANCE relative to oxygen 
is determined. Because of 
the contribution of nitrogen 
to the oxygen track, which is 
in the order of 10$, the 
abundance ratios are corrected 
by 10 % compared to the pure 
number ratios. The result 
is shown in table 1 along 
with experimental results 
of other authors. The experi- 
mental errors are based on the 
statistical errors and en- 
larged by an estimated upper 
limit of systematic error 
due to uncertainties in the 
energy calibration. 



’ KINETIC ENERGY (MeV/NUCL.) 


The ENERGY DEPENDENCE 
cf the abundance ratios for 
carbon and iron have been 
determined. Figure 4 shows 
the result along with 
similar ratios obtained by 
other experimenters. 


Fig. 3. Energy spectra of 
carbon, oxygen and the iron 
group derived from figure 2. 
To get the fluence of the 
flare multiply by 60. 


Our data show a carbon/oxygen ratio of about 0.4, which is 
nearly independent of energy and which agrees reasonably well 
with the value of other authors at higher energies. The results 
on iron show an enhanced flux compared to the results at higher 
energies of Teegarden et al . , and Bertsch et al . . We do not see 
any significant change in the ratio over the energy range from 
0.62 to 6.9 MeV/nucleon, which seems to be in contradiction to 
Braddy et al. . 1973, which is the only other data; overlapping 
in energy with ours. 
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Table 1 Relative Abundance ratios (Oxygen - 1.0) 



Energy 
MeV/Nucl . 

Present Work 

Crawford 
et al. 1972 

Tocgarden 
et al.1973 

— 

Bertsch 
et al.1972 

Mogro-Con- 
ero et al.1972 

m 

0.62-4.21 

0.43+0.02 


0.49+0.03^ 

0.56+0.06 

- 


0.62-6.9 

1.00 

1.0 

1.0 

1.0 

1.0 

Ne 

0.73-6.9 

0.13+0.02 

_ , + .l 
0,1 -.05 

.127+.011* 2 ^ 

.16+. 03 

.21+. 08 

Mg 

0.96-6.9 

0.22+0.02 

. 125+.03 

.182+.014 (2) 

.056+. 014 

.13+. 06 

SI 

0.89-6.9 

0.23+0.03 

.1 

.107+. Oil * 2) 

.028+. Ol 

.36+. 12 

Fe 

0.62-5.2 

0.17+0.02 

.09 +.03 

.028+.005* 2 * 

•06+.015 <33 

.37+. 12 


( 1 ) Original data was normalized to SI - l.(0/Si - 10.0). (2) Sept. 1.1971. (3) Bertsch et al.1973 


Pig. C/0 and Fe/O 

abundance ratio as a 
function of energy. The 
source of the data are 
indicated in the figure. 
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Discussion. From 

table 1 we see a reason- 
able good agreement of 
the C/0, Ne/0 and Mg/0 
ratios with the values 
of Teegarden e t al. , 

1973. The data of 
Bertsch et al. , 1972 
as measured at consider- 
ably higher energies 

E> 20 MeV/nucleon are a factor of five for Mg/0 and a factor of 
eight for Si/0 lower than our values. 
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The values of Crawford for Mg, Si and Fe are about a 
factor of two lower than ours. This difference may in part be 
explained because their data originally was normalized to 
silicon and had to be converted to the value relative to oxygen, 
which introduces additional uncertainties in the absolute 
values. On the other hand their data overlap in energy only at 
the lowest energy point (3 MeV/nucleon). Very recently Braddy et 
al. , 1973 published data of the Fe/O ratio which completely over- 
lap our energy range. Their Fe/O ratio is about a factor 2 lower 
than ours and shows a fairly strong energy dependence in a 
region, where we see no such differences in the abundance ratios. 
It should be emphasized that their and our measurements are not 
affected by the magnetosphere, contrary to rocket measurements 
at these low energies. 

Analysing several solar flare events recently Mogro-Com- 
pero et al. , 1972b and 1973b found. a variability in the Fe 
abundance ratio which is correlated to the CNO/He ratio, in the 
sense that a He rich event also is rich in iron relative to the 
C, N, 0 nuclei. The authors relate the variation to partial 
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ionization during the acceleration process. If this is true and 
provided that the charge equilibrium is maintained during the 
acceleration process, one would expect the strongest energy de- 
pendence of the ratio in the energy region, where the ionic chaige 
to mass ratio Q/A of iron and oxygen changes most. But for the 
October 72 event we see no energy dependence of the Fe/O ratio 
Just in this energy region between few hundred keV/nucleon and 
several MeV/nucleon. 

Our results would suggest just the opposite that at least 
for this flare during the acceleration process Q/A is the same 
for both elements and independent of energy. A possibility would 
be that the particles get accelerated after being fully stripped 
by energetic electrons in hot flare regions, as proposed by 
Cartwright et al. , 1972 for the 'first step' of their accelera- 
tion mechanism. This would imply that the acceleration takes 
place in a region of the solar atmosphere high enough that the 
density is sufficiently low that no charge exchange by collision 
alters the state of ionization and the energy of the particles 
before escape into the interplanetary space. 
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